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1 INTRODUCTION 

The belief that molecular clouds are gravitationally bound 
objects has led to a long-standing problem in star formation. 
If clouds are bound, and if they collapse in of order a free-fall 
time, then the rate of star formation should be around two 
orders of magnitude greater than what is observed 



|erman fc Evans||1974| . To circumvent this there have been 
two main approaches. Firstly, one can assume that molecular 
clouds undergo collapse on longer timescales, requiring that 
some process prevents, or at least slows their collapse. This 
could be magnetic fields and slow ambipolar diffusion ( |Shu| 
et al.|[l987l |Basu Mouschovias|[T994l | Allen Shu||2000| 
Mouschovias et al. 2006; Sh u et al.|2007|, or mic roturbulence 



ABSTRACT 

The most recent observational evidence seems to indicate that giant molecular clouds 
are predominantly gravitationally unbound objects. In this paper we show that this is 
a natural consequence of a scenario in which cloud-cloud collisions and stellar feedback 
regulate the internal velocity dispersion of the gas, and so prevent global gravitational 
forces from becoming dominant. Thus, while the molecular gas is for the most part 
gravitationally unbound, local regions within the denser parts of the gas (within the 
clouds) do become bound and are able to form stars. We find that the observations, in 
terms of distributions of virial parameters and cloud structures, can be well modelled 
provided that the star formation efficiency in these bound regions is of order 5-10 
percent. We also find that in this picture the constituent gas of individual molecular 
clouds changes over relatively short time scales, typically a few Myr. 

(e.g. Bertoldi fe McKee"||1992| |Dib et al.||2007| ), where a v is 
the line of sight velocity dispersion and R is a measure of the 
radius of the cloud. If the cloud is in virial equilibrium, then 
a — 1 and 2T+VK=0 where T is the kinetic and W the grav- 
itational energy, whereas T + W = corresponds to the zero 
energy configuration. The value of a is simply a measure of 
the ratio of kinetic to gravitational energies, and finding a 
both observationally and from simulations is highly uncer- 
tain, depending on the mass and radius determinations, and 
typically does not account for magnetic fields, surface terms 
or pr ojection effects ([B allesteros-Paredes 2006] |Dib et ah] 
2007| |Shetty et al.|2010| ). Thus the estimated value of a for 
one cloud may not predicate its consequent evolution, but 



(Zuck- 



(Krumholz & McKee 



2005 



Krumholz et al. 



2006). How- 



ever turbulence can 


also promote collapse on small scales 


(Klessen et al. 2000 


|Mac Low & Klessen||2004|), (includ- 


ing MHD turbulence 


Heitsch et al.|2001||Vazquez-Semadeni 


et al.|2005 Elmegreen|2007|. Alternatively, in light of recent 



the distribution of a gives an indication of the importance 
of gravity in a population of clouds, and whether they are 
predominantly bound or unbound. 



observational evidence, one can assume that clouds are typ- 
ically short -lived entities (|Hartmann et al.|2001 |Elmegreen| 
2002[|2007||Ballesteros-Paredes et al.|2007| ) and posit mech- 



In a recent observational study, Heyer et al. (2009) pub- 
lished revised estimates of molecular cloud masses, sizes 
and virial parameters from the previous seminal work by 



Solomon et al.| ( |1987| ) . Although they claim that molecular 



anisms which prevent most of the gas forming stars, such 



clouds are virialised, their plots seem to indicate that in fact 



as stellar feedback and magnetic fields ( Vazquez- Semadeni 
et al.||2005| |Elmegreen||2007| |Price fe Bate||2008) . However 



most of the clouds are unbound (as also found by Heyer et al 
2001 ). In Fig. 1 (lower middle panel) we plot the value of the 



if the cloud is not globally bound in the first place, as sug- 
gested by recent observational evidence, we already amelio- 
rate the disparity between the existence of short-lived molec- 
ular clouds and the global star formation rate. 

The virial parameter of a molecular cloud is usually 
defined as 



virial parameter taken from the clouds observed by |Heyer| 
et al. (2009). Heyer et al.| ( |20C)9] ) suggest that the masses 



5ajR 
GM 



(1) 



of their clouds are underestimated by a factor of 2 due to 
non-LTE effects and CO abundance variations (see also cal- 
culations by |Glover et ah] ( |2010| ) and |Shetty et al] (2010) ) 
within a cloud, so we have doubled their cloud masses to 
produce the panel in Fig. 1 (lower middle panel). It is evi- 
dent that most of the observed clouds have virial parameters 
larger than unity, indicating that most clouds are not gravi- 
tationally bound. Even if we take a > 2, whereby the clouds 



2 C. L. Dobbs, A. Burkert and J. E. Pringle 



are strictly unbound, this still leaves 50 per cent unbound 
clouds. Similar distributions of a are also found for exter- 
nal GMCs, including those in M31 ( |Rosolowsky||2007[ from 
a sample of 10 5-6 M clouds), and the GMCs detected in 



2 SIMULATIONS 



several nearby galaxies by Bolatto et al. ( 2008 ) 



In a recent paper Balle steros-Paredes et al.| (|2010) sug- 
gest that GMCs are undergoing hierarchical gravitational 
collapse, whereby the collapse occurs on scales from indi- 
vidual cores to the whole cloud. However it is not neces- 
sary that the cloud should be globally gravitationally bound. 
Simulations of unbound, turbulent, GMCs naturally lead to 
localised star formation, rather than spread over the entire 
cloud (|Clark et al.||2005| |2008|. This then naturally leads 



to a lo w star formation efficiency. Recently, |Bonnell et al.| 
(2010) performed calculations of an unbound 10 Mq cloud, 



and showed that stellar clusters form in bound regions of the 
cloud. The internal kinematics of these clouds could be due 



to cloud-cloud collisions or large scale flows (Bonnell et al. 



2006l|Dobbs fe Bonnell|[2007] |Klessen fe Hennebelle||2010| ), 

and/or stellar feedback (e.g. |Mac Low &; Kle ssen 2004 and 
references therein). 

Most numerical work has tended to focus on calculat- 
ing the virial parameters of clumps within giant molecular 
clouds ( |Dib et al.|2007||Shetty et al.|2010[ ), which, since they 
are the sites of star formation in GMCs, are more likely to 
be bound. However in recent simulations of a galactic disk, 
it has been possible to identify individual GMCs and de- 
termine their virial parameters (Dobbs 2008; Tasker & Tan 
2009). These results show that the virial parameter, a (see 



Section 2.2) typically lies in the range of around 0.2 to 10. 
In this paper, we address the question of how molecular 



clouds can remain unbound. Pringle et al. ( 2001 ) argued that 
if molecular clouds are short-lived, with lifetimes comparable 
to a few tens of Myr, then they must be formed from a large 
reservoir of dense interstellar gas, which may or may not 



itself be molecular. Dobbs et al. (2006) has shown that the 



formation of the global structure of molecular gas (clouds, 
spurs etc.) does not in itself require self-gravity, but that 
formation can come about for entirely kinematic reasons. In 
this paper we take these ideas a step further and attempt to 
model the observed properties of molecular clouds. We self 
consistently follow the evolution of clouds in a galactic disc, 
taking into account cloud collisions and cloud dispersal by 
energy input from stellar feedback. The clouds we consider 
are of size tens of parsecs, we are unable to resolve very small 
clouds. The particular properties we try to match are the ob- 
served distribution of the virial parameter a, the shapes of 
the clouds and their internal structures. We find that these 
properties can be matched simply by assuming that those 
regions within molecular clouds that become self- gravitating 
are able to form stars at some small efficiency (5-10 per 
cent) which gives rise to feedback in the form of input of en- 
ergy and momentum (Section 2). Thus if say only around 10 
per cent of a cloud is bound at any one time, and those parts 
form stars at around 10 per cent efficiency, the problem of 
the two order of magnitude difference in the star formation 
rate identified by |Zuckerman &; Evans (1974]) can be over- 
come (see Section 3). We demonstrate that with this simple 
assumption, those structures which would be identified as 
molecular clouds are, for the most part, globally unbound, 
with properties giving a reasonable match to the data. 



The calculations presented here are 3D SPH simulations us- 
ing an SPH code developed by Benz ( Benz et al.|1990 ), Bate 
( |Bate et al.|1995| ) and Price ( |Price fe Monaghan|2007| . The 
code uses a variable smoothing length, such that the density 
p and smoothing length h are solved iteratively according to 



Price & Monaghan ( 2007 ) , and the typical number of neigh- 
bours for a particle is ~ 60. Artificial viscosity is included 
to treat shocks, with the standard values a = 1 and /3 = 2 
(Monagh an|1997 ). In all the calculations presented here, the 
gas is assumed to orbit in a fixed galactic gravitational po- 



tential. The potential includes a halo (Caldwell & Ostriker 
|1981| ), disc ([Binney &; Tremaine||1987 ) and 4 armed spiral 



component ( |Cox &; G6mez|2002 ). The gas particles are ini- 
tially set up with a random distribution, and assigned veloc- 
ities according to the rotation curve of the galactic potential 
with an additional 6 km s _1 velocity dispersion. The rota- 
tion curve is flat across most of the disc, with a maximum 
velocity of 220 km s _1 . 

We present results from 4 different calculations, as sum- 
marised in Table 1. Run A was already presented in [Dobbs 
( |2008| ), and is more simplistic than Runs B, C and D. The 
total gas mass is 5 x 10 9 M in Run A, and 2.5 x 10 9 M in 
Runs B, C and D, and corresponds to one or two per cent 
of the total mass of the galaxy. The surface densit y of the 
Milky Wa y is about 12 M pc -2 (it is 10 M pc -2 in |Womre| 



et al.|2003 excluding helium) thus a little higher than Runs 
B, C and D. The mass resolution is 1250 M for Run A, 
and 2500 M for Runs B, C and D. 

For Run A, we allocate particles at radii between 5 and 
10 kpc. The gas is assumed to be a two phase fluid. The in- 
terstellar medium has two isothermal components, one cool 
and one warm. We omit thermal considerations and so there 
is no transition between the two phases; the cool gas remains 
cool and the warm gas remains warm, throughout. The cool 
and warm gas comprise equal mass in the simulations. The 
gas has initial scale heights of 150 and 400 pc in the cold and 
warm components respectively, but these decrease to 20-100 
pc and 300 pc with time. The mean smoothing length is 40 
pc. In Run A, we also include a magnetic field, such that 
the plasma /3 of the cold gas is 4. The magnetic field is im- 
plemented using Euler potentials as described in Dob bs <£] 



Price (2008) 



In the remaining calculations (B, C and D), we inves- 
tigate the effect of stellar feedback. In these cases we allow 
the ISM to exhibit a multiphase nature from 20 K to 2 x 10 6 
K. The cooling and heating of the ISM is calculated as de- 
scribed in |Dobbs et al.| ( |2008| ) . Apart from feedback from 
star formation, heating is mainly due to background FUV, 
whilst cooling is due to a variety of processes including colli- 
sional cooling, gas-grain energy transfer and recombination 
on grain surfaces. The gas initially lies within a radius of 
10 kpc, and has an initial temperature of 7000 K. The im- 
plementation of stellar feedback will be described in detail 
in a forthcoming paper, but a simple description is included 
here. The gas is assumed to form stars when a number of 
conditions are met, i) the density is greater than 250 cm -3 , 
ii) the gas flow is converging, iii) the gas is gravitationally 
bound (within a size of about 20 pc, or 3 smoothing lengths), 
iv) the sum of the ratio of thermal and rotation energies to 
the gravitational energy is less than 1, and v) the total en- 
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ergy of the particles is negative (see Bate et al. 1995). If 



all these conditions are met, we assume that star formation 
takes place; there is no probabilistic element in our calcula- 
tion. We do not however include sink particles, instead we 
deposit energy in the constituent particles. We present cal- 
culations with star formation efficiencies, e, of 1, 5 and 10 
per cent (Runs B, C and D respectively). This means that of 
the mass that satisfies the above criteria, a fraction e of the 
molecular gas contained therein is assumed to form stars 
instantaneously and to provide an energy input (approxi- 
mately 1/3 thermal and 2/3 kinetic energy) of 10 51 ergs per 
160 M of stars assumed to form.^This energy input, com- 
bined with our cooling and heating prescription, leads to a 
multiphase ISM. In the case of Run C (e = 5 per cent), from 
150 Myr onwards approximately one third of the gas lies in 
the cold, unstable and warm regimes. 



2.1 Locating clouds 

We identify clouds using the same method as described in 
Dobbs (2008). We apply a clumpfinding algorithm, which 
simply divides the simulation into a grid, and locates cells 



over a given surface density. Adjacent cells which exceed 
this criterion are grouped together and labelled as a cloud. 
Clouds which contain less than 30 particles are discarded, 
thus clouds in Runs B, C and D are at least 7.5 x 10 4 Mq 
(and clouds in Run A 3.75 x 10 4 M ). The mean number 
of particles in a cloud is ~ 85 for Runs A, B and C. The 
properties of the clumps reflect the total, rather than the 
molecular gas, but we would typically expect these clumps 
to exhibit high molecular fractions. For most of the re- 
sults we present, we chose a surface density criterion of 100 
M©pc -2 ~ 4 x 10 22 cm -3 . Changing this criterion has little 
effect on the distribution of the virial parameter, it merely 
reduces or increases the number of clouds selected. 

For Runs C and D (e = 5, 10 per cent), we are able to 
run the simulation for sufficiently long (300 Myr) that we 
can calculate cloud properties when the system has roughly 
reached equilibrium. However for Runs A (high surface den- 
sity, no feedback) and B (1 % efficiency) we are limited by 
the high surface densities reached by a large fraction of the 
gas. 

To illustrate the global structure of the disc in our mod- 
els, the column density of the gas in Run C (e = 5 per cent) 
is shown at a time of 200 Myr in Fig. 2. The dense gas is ar- 
ranged into clouds along the spiral arms and spurs extending 
from the arm to interarm regions. 



2.2 Virial parameter 

In determining the virial parameters for our clouds, we cal- 
culate a as shown in Eqn. 1, where a v is the line of sight 
velocity dispersion and R is defined as the radius of a circle 
with the equivalent area of the cloud. This corresponds to 



that used by Heyer et al. (2009). We take bound clouds as 
having a < 1. 

In the case with magnetic fields (Run A) , we find largely 



1 This corresponds to a Salpeter IMF with limits of 0.1 and 100 

M©. 
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Figure 2. The gas column density is shown for Run C (e = 5 per 
cent) at a time of 200 Myr. Dense gas, corresponding to the clouds 
located in the analysis presented here, predominantly lies along 
the arms, and spurs which extend from the arms into interarm 
regions. 



unbound clouds, where local gravitational collapse is pre- 
vented by magnetic pressure. With minimal stellar feedback 
(Run B, e = 1 per cent), we obtain many more bound clouds, 
particularly at higher masses. This clearly disagrees with 
the observations. For both Runs C and D (e = 5 and 10 
% respectively), we find that the clouds are predominantly 
unbound, and the distributions in the virial parameter, a 
are in agreement with the observations. This can be seen 
visually and is confirmed by comparing the distributions of 
a using the KS test (see also Fig. 1). Given the uncertainties 
in determining a, if we require a > 2 for an unbound cloud, 
about half of the clouds in Runs C and D are unbound. 
There is little change in the fraction of unbound clouds with 
mass (and therefore resolution) in these calculations, with 
the exception of Run D, where feedback is responsible for 
preventing bound, higher mass clouds. 



2.3 Evolution of individual clouds 

The evolution of an individual cloud is often very complex, 
involving collisions, fragmentation and dispersion by feed- 
back. Moreover the gas which constitutes a cloud can change 
on relatively short timescales. Thus studying the evolution 
of individual clouds, and establishing why they never, or 
rarely become gravitationally bound is not straightforward. 
In this Section we illustrate this behaviour by studying the 
nature and development of individual GMCs in the different 
calculations. 



2.3.1 Cloud development in Run A 

In Fig. 3, we highlight the contribution of collisions to the 
internal velocity dispersions of clouds in Run A. We show a 
collision between two clouds in Run A (which includes mag- 
netic fields), where small scale gravitational collapse does 
not occur. After the collision between the two clouds, the 
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Run Surface density 


ISM 


P 


e 


No. particles 


Time chosen to 


(M pc- 2 ) 






per cent 




locate clouds (Myr) 


A 20 


Two phase isothermal 


4 


N/A 


4 x 10 6 


130 


B 8 


Multiphase (20 - 2 x 10 6 K) 


oo 


1 


10 6 


110 


C 8 


Multiphase (20 - 2 x 10 6 K) 


oo 


5 


10 6 


200 


D 8 


Multiphase (20 - 2 x 10 6 K) 


oo 


10 


10 6 


200 



Table 1. The different calculations performed are described above. Run A was presented in Dobbs (2008). For this two phase isothermal 
calculation, half the gas lies in the warm (10 4 K) phase whilst half is cold (100 K). e is the star formation efficiency in the calculations with 
feedback (see text for details). For Runs C and D, the time we locate the clouds is not important as they have reached an approximate 
equilibrium state. 





Figure 1. The distribution of the virial parameter (a) is plotted with mass for clouds identified in Run A (top left, with magnetic 
fields), in the calculations with feedback adopting efficiencies of 1, 5 and 10 %, and the |Heyer et al.| ( |2009| data (lower middle). We find 
a population of predominantly unbound clouds, in rough agreement with the observations, for the models where localised gravitational 
collapse is limited by magnetic fields (Run A), or gravitational collapse occurs but there is a realistic level of stellar feedback (Run C, 
top right, Run D, lower left). There are many more bound clouds for the case with a very low level of stellar feedback (Run B, top 
middle). In the final panel (lower right), the cumulative fraction of clouds with a given a is shown for the Heyer data (dotted) and for 
the clouds from Run C, with 5 per cent efficiency (solid line). The KS test confirms that the distributions of a from the observations 
and simulations match, giving P = 0.11 and P = 0.21 for Runs C and D respectively. 



velocity dispersion increases, which means the virial param- 
eter also increases. The increase in the velocity dispersion is 
prolonged because the clouds contain substantial substruc- 
ture - the merging of this substructure is seen in the mid- 
dle panel. Thus the collisions between clouds are not really 



dissipative (as stated in Dobbs et al. (2006)); rather the 
energy is transferred to the internal motions of the clouds. 
We also simulated the cloud interaction shown in Fig. 3 in 
isolation, and at higher resolution, without magnetic fields 
or feedback. This confirmed that the collision induces ran- 
dom large scale motions which prevent widespread collapse 



in the cloud for around 10 Myr, independent of the mag- 
netic field. This demonstrates that the energy input from 
cloud-cloud collisions can be comparable to, or even exceed 
the energy dissipated, i.e. a collision can lead to an increase 
in the global virial parameter. The generation of random 
velocities is analogous to that presented in previous work 
dBonnell et aLl|2006| [Dobbs fe Bonnell||2007| ), and relies on 
the assumption that the ISM is clumpy on all scales. 
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Figure 3. The evolution of a cloud from Run A is shown at 
times of 123 (top), 132 (second) and 141 (third panel) Myr. The 
evolution of the mass (in units of 10 6 Mq), velocity dispersion and 
virial parameter are shown on the lowest panel. The two clouds 
in the top panel merge, and there is a subsequent increase in both 
the velocity dispersion and a. They form a single clump of mass 
7 x 10 5 Mq. The initial clumps are coloured blue and red, so the 
constituent particles can be traced at later times. The particles 
coloured green are particles present in the clouds at 132 and 141 
Myr, which were not present in the clumps at 123 Myr. 
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Figure 4. The evolution of a cloud from Run C (with 5 % effi- 
ciency stellar feedback) is shown, at 189 (top), 198 (second), 200 
(third) and 202 (fourth panel) Myr. The cloud is formed by the 
merger of smaller clumps. Stellar feedback events (for example the 
cross in the second panel) then alter the shape of the cloud and 
finally result in the separation of the cloud into several separate 
clumps. Separate clumps, (as picked out by the clumpfmding algo- 
rithm), are shown simply in different colours, but the constituent 
particles are not all the same at different times. For example only 
2/3 of the particles in the cloud at 198 Myr are in the cloud shown 
at 200 Myr. The right hand panels show column density images 
of the clumps and their surroundings. The white boxes indicate 
the size of the regions shown on the left hand panels. Fig. 5 shows 
the evolution of a, the mass and the velocity disperions, a, of the 
cloud, and the constituent clumps which formed the cloud. 
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Figure 5. The evolution of a (green, solid), mass (red, dashed 
in units of 10 6 Mq) and a (blue, dotted) of the clouds shown 
in Fig. 4. Cloud-cloud interactions have some role in determining 
the dynamics of the cloud, and maintaining cr, but a is dominated 
by stellar feedback. The different lines at 189 Myr (and the line 
at 193 Myr) correspond to several clumps merging to form one 
large cloud (193 Myr to 200 Myr). After 200 Myr, the cloud again 
splits up into multiple components, due to feedback. 



2.3.2 Run C: a multiple cloud interaction with feedback 

In Fig. 4 we show the evolution and interaction of a multiple 
set of clouds in Run C (e =5 per cent). A single cloud was 
selected at a time of 198 Myr, and the clouds which contain 
the same particles were identified at earlier and later times. 
We find that the cloud identified at 198 Myr is formed from 
the mergers of several smaller clouds. In the first panel (189 
Myr) we can identify 5 separate clouds. As these merge to 
produce a cloud of 2 x 10 6 Mq some 9 Myr later, the effects of 
stellar feedback can be seen for example in the cloud in the 
second panel (a bubble blown out by feedback is indicated 
by the cross). Feedback plays a large role in shaping the 
cloud, and regulating the dynamics. The clouds, as picked 
out by the clumpfinding algorithm (left hand plots) show 
much more filamentary structures compared to the clouds 
in Run A (Fig. 3). By 202 Myr (4th panel), stellar feed- 
back has succeeded blowing away the top part of the cloud, 
and splitting the cloud apart. Over the course of the plots 
shown (13 Myr), there are 5 supernovae events in the cloud. 
In Fig. 5 we show the corresponding evolution of a and the 
velocity dispersion. It can be seen that the velocity disper- 
sion is maintained at about 6 km s _1 , and the virial parame- 
ter, a, above unity throughout. Thus in this case, a multiple 
cloud collision together with feedback (from small regions 
of the cloud which become bound and allow star formation) 
maintains the unbound nature of the cloud as a whole. 

It can also be seen from Fig. 4 (right hand panels) that 
the clouds we identify are part of a larger region of dense 
gas, which is hundreds as opposed to tens of parsecs in size. 
In a galaxy with a high molecular gas fraction such a feature 
would correspond to a Giant Molecular Association (GMA). 
Whilst these regions are still unbound in our calculations, 
they appear to have a longer lifetime than the GMC sized 
clouds. 
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Figure 6. The evolution of a second cloud from Run C, with 
5 % efficiency stellar feedback. Unlike the cloud shown in Fig. 4, 
this cloud is too massive to be disrupted by feedback or collisions 
and becomes increasingly more massive, and bound with time. At 
160 Myr (top panel), the cloud is still filamentary, and marginally 
unbound. By 200 Myr (lower panel), the cloud contains no fila- 
mentary structure and is strongly bound. The right hand panels 
show column density images of the regions containing the clouds 
on the left, the white boxes indicating the sizes of the left hand 
plots. The lowest panel shows the evolution of en, mass (in units 
of 10 6 M ) and a for the cloud. 



2.3.3 Run C: the evolution of an isolated massive cloud 

For the calculation, with 5 % efficiency feedback (Run C), 
the timescales for the majority of clouds to merge and be- 
come disrupted are relatively short, of order several Myr. 
The exceptions are two longlived bound clouds, which have 
masses of 3 x 10 6 M and 5 x 10 6 M respectively. These are 
the high mass points seen in Fig. 1 (top right plot). We show 
the evolution of the 5 x 10 6 Mq cloud in Fig. 6 over a period 
of 40 Myr. The top panel shows the cloud at a time of 160 
Myr, when the cloud is clearly irregular in shape. By 200 
Myr, the cloud has a much more regular, quasi-spherical ap- 
pearance and is not filamentary in any way. This cloud finds 



Why are most molecular clouds not gravitationally bound? 7 



itself in between the spiral arms, and does not undergo any 
significant interactions with other clouds after 160 Myr. In 
the lower panel, we plot the evolution of a, the velocity dis- 
persion and the mass. We see that the cloud is continuing to 
accrete material, and grow in mass, becoming steadily more 
bound. Feedback (with e =5 per cent) is insufficient to dis- 
rupt the cloud, although feedback does maintain a constant 
velocity dispersion of ^6 km s _1 . It is likely that this cloud 
would eventually form a bound stellar cluster, though we 
do not attempt to follow this in our calculation. In Run B 
(e = 1 per cent) there are many more clouds which display 
this behaviour. 



10 b 

Mass (M ) 



2.4 The constituent gas of the clouds 

In the current paradigm of molecular cloud formation and 
evolution, GMCs are assumed to be bound objects which 
consist of essentially the same gas for the duration of their 
lifetimes. In Fig. 7 we take all the clouds at a given time 
in Run C (e = 5 per cent), and plot the percentage of gas 
which remains in a given cloud after 10 Myr. 

In all cases we find that the constituent gas of the clouds 
changes on timescales of < 10 Myr. We find that about 
50 per cent of clouds are completely dispersed within 10 
Myr. A substantial fraction of clouds are shortlived, either 
dispersing to lower densities, merging with other clouds to 
produce more massive clouds, or some combination of these 
processes. There are a few clouds which substantially retain 
their identity over a period of >10 Myr. 

This highlights that generally the constituent gas in 
GMCs is likely to change on timescales of Myr. This may 
mean that discussing clouds lifetimes, which are thought to 
be 20-30 Myr ( |Leisawitz et al.|1989||Kawamura et al.|2009| ), 
may not make sense. A cloud seen after 30 Myr may not be a 
counterpart to any cloud present at the current time. In our 
calculations this is only true for the most massive clouds. 

Thus we see that the clouds tend to display a variety 
of behaviours. The relatively low mass GMCs undergo fre- 
quent collisions, are readily disrupted, and a will change 
accordingly. Even if the cloud becomes bound, it may un- 
dergo another dynamical interaction on a short timescale, 
and become unbound. In contrast the more massive clouds 
undergo less dramatic behaviour. There are relatively few 
clouds of this mass, so they very rarely undergo collisions 
with objects of a similar size, whilst above a certain mass 
they are not so easily torn apart by feedback. 



2.5 The shapes of clouds 

From our calculations with different levels of feedback, we 
obtain distributions of clouds which are predominantly un- 
bound (Runs C and D, with 5 and 10 % efficiency) or bound 
(Run B, with 1 % efficiency). We have demonstrated that the 
observations are most likely fit by a distribution of mainly 
unbound clouds. 

In addition, we find that, the bound clouds are much 
more regular, spherically shaped, whilst the unbound clouds 



Figure 7. The percentage of gas which still lies in the same 
cloud after 10 Myr is plotted versus the cloud's mass for clouds 
identified in Run C (with stellar feedback and e = 5 per cent). 
This percentage is calculated by locating clouds at 2 timeframes, 
10 Myr apart. We search for the constituent gas particles of a 
given cloud in the clouds present 10 Myr later. Sometimes there 
may be more than one cloud containing the particles of an earlier 
cloud (and in some instances no clouds!), in which case we select 
the cloud which has the maximum of the particles the same as 
the cloud at the earlier time frame. In 35 per cent of cases the 
clouds are completely disrupted, whilst the median amount of gas 
remaining in the cloud is 22 per cent. 



1.5 and 2. We also show in Fig. 8 the distributions of aspect 
ratios for the calculations with low feedback Run B (pre- 
dominantly bound clouds), and the higher feedback case, 
Run C (predominantly unbound clouds). 

In Run C (centre panel), with feedback efficiency of 5 
per cent, the distribution of aspect ratios has reached an 
equilibrium state, and is found to be similar to observations, 
with a peak at about 1.5. In Run B (left panel), with 1 % 
efficiency, equilibrium has yet to be achieved, and more and 
more clouds have aspect ratios of around unity with increas- 
ing time. As we would expect, the virialised clouds tend to 
have aspect ratios close to unity. The prominent peak at as- 
pect ratios of unity does not agree with the observations, 
reconfirming our previous conclusions that the simulations 
with efficiencies of 5 - 10 per cent are in best agreement 
with the observations. 



3 DISCUSSION 

In this paper we have addressed the recent observational 
evidence that most molecular clouds within the Galaxy are 
not gravitationally bound. This evidence contrasts with the 
original claims of 



Solomon et al. 



(1987), long since propa- 



have very irregular shapes. Koda et al. (2006) carried out a 



study of Galactic molecular clouds and determined the as- 
pect ratios of the clouds. Their results are shown in Fig. 8, 
and indicate the most common aspect ratios are between 



gated into the field, that molecular clouds are bound and in 
virial equilibrium. The idea that molecular clouds are bound 
entities has also been taken as a starting point for many the- 
ories of star formation. 

Of course, for star formation to take place it is neces- 
sary that some parts of a cloud be self- gravitating and able 
to undergo collapse. What the observations seem to indicate 
however is that it is not necessary for the cloud as a whole to 
be dominated by gravity. With this in mind, we present here 
simulations of the ISM in a galaxy with a fixed spiral poten- 
tial. By including simple heating and cooling of the gas, we 
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Figure 8. The distribution of aspects ratios of the clouds is shown for the models with 1 % efficiency feedback (Run B, left), and 5 
% efficiency feedback (Run C, centre). The distribution of aspect ratios for Galactic clouds is shown on the right ( |Koda et al.||2006) . 
The clouds for the 5 % efficiency case (centre) reasonably match the observations, although even in this case our clouds are slightly 
more peaked towards low aspect ratios than the observations. The distribution does not change with time, once equilibrium has been 
established. With 1 % efficiency (left), the distribution evolves to a strong peak at 1, in definite contradiction to the observations. 



are able to identify those parts of the ISM which are dense 
enough to represent molecular gas, and so are able to iden- 
tify what would be observed as molecular clouds. We allow 
the parts of those clouds which are sufficiently dense and 
sufficiently gravitationally bound to notionally form stars. 
Because such clouds are generally highly inhomogeneous en- 
tities, within them there will be some regions (in our simu- 
lations typically representing only < 30 percent of the mass) 
which are gravitationally bound, and within which star for- 
mation takes place. This star formation is taken to manifest 
itself as an input of energy and momentum into the sur- 
rounding gas. The global galactic star formation rate is in 
accordance with the results of |Kennicutt| ( [T998| ). 

Using this simple, and highly idealised, input physics 
we are able to reproduce both the observed distribution of 
virial parameters of molecular clouds in the Galaxy (with 
most of the clouds being unbound) and also the observed 
distribution of cloud shapes (in terms of their aspect ra- 
tios). We find that the velocity dispersions within clouds 
are maintained not just by feedback from star formation 
but also by collisions between non-homogeneous clouds (cf. 
Dobbs fe Bonnell|2007| ). However with no, or little feedback, 
the clouds are predominantly bound and quasi-spherical (as 
found in Run B and by Tas ker &; Tan|2Q09 ) , in disagreement 
with observations. 

We also find that the constituent material of a typi- 
cal cloud only remains within that cloud for a timescale of 
around 10 Myr. Thus for timescale much longer than this, 
the concept of a cloud lifetime is no longer meaningful. 

We note that the properties and lifetimes of clouds de- 
pend somewhat on the size scales considered. Above some 
surface density threshold, we would expect to start select- 
ing bound regions within a GMC, and therefore we would 
obtain a higher fraction of bound objects. We have not con- 
sidered the properties of larger GMAs either. The fraction 
of unbound clouds also depends on how we define a, and 
what threshold we use. However we note that for our clouds, 
even if a is low, external energy input from collisions, and 
feedback sources within a cloud can act to increase a v , and 
therefore a. Thus the main difference between our mod- 
els and previous analaysis, for example that presented by 



Ballesteros-Paredes et al. (2010), is that collisions and feed- 



back play a much more important role. 

In summary, the idea that all molecular clouds are grav- 
itationally bound entities is neither observationally viable, 
nor theoretically necessary. It is no real surprise that most 
molecular clouds identified in the Galaxy are globally un- 
bound, and that the rest are at most only marginally bound. 
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